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gluon field fluctuations, within the context of the stochastic vacuum 
e-j ! model (SVM) of pure QCD. After formulating the problem in terms 

of the generating functional of field strength cumulants, we argue that 
the resulting predictions about the Casimir energy are compatible with 
the phenomenologically required bag energy term. 
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Abstract 



1 Introduction 

In order to adjust hadron spectroscopy, the MIT bag model requires a bag 
energy term of the form —Z/R with a positive Z of order one [T]. One of 
the main contributions to this energy is naturally ascribed to the Casimir 
energy of the fields involved. However, due to the daunting difficulty posed 
by Casimir energy calculations in nonabelian gauge field theories, this prob- 
lem has only been approached under rather strong simplifying assumptions. 
Indeed, vacuum effects coming from the pure Yang-Mills sector have been 
computed by considering a simplified description, which treats gluon modes 
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as a collection of free, photon-like gauge fields. This assumption (usually 
justified by invoking asymptotic freedom) leads, however, to an energy with 
the wrong sign: ~ +0.7 /R. Thus the effect becomes repulsive, instead of the 
attractive one needed to comply with phenomenology [2J-(S|- For a discussion 
of the bag model, as well as of the current status of Casimir calculations in 
the bag, see, for example, refs. |7| and [8]. 

Neither the inclusion of the (rather small) fermionic Casimir energy con- 
tribution ~ — 0.02/ R due to three light quark modes [9], nor the (relevant) 
effect of the quarks' center of mass quantum fluctuations, which is attractive, 
~ —Z cm /R, Z cm ~ 0.6 — 0.8, are sufficient to reach the required value of Z 
(see the discussion in [7] and [TO]). 

It has been pointed out that a resolution to this problem could be associ- 
ated with the intrinsically nonabelian nature of the Yang-Mills theories, and 
their ensuing nonperturbative effects. This is certainly plausible, although it 
is not at all evident how to include those effects in a tractable way. 

In this article, we present an approach that enables one to take partly 
into account the nonperturbative dynamics of QCD into the calculation of 
the Casimir energy, and argue that its predictions are phenomenologically 
sound. 

We first show how the Casimir energy in pure QCD can be written in 
terms of a subclass of gauge invariant correlators, namely, the shifted field- 
strength correlators. Lattice simulations show [11] that these correlators 
exhibit Gaussian dominance, so that they can be described in terms of the 
so called stochastic vacuum model (SVM) of pure QCD [12] (for a review 
see refs. [33] and [H]). This model can be parametrized in terms of just two 
structure functions, obtained by fitting the lattice data, and it consistently 
describes short distance perturbative aspects, such as asymptotic freedom, 
as well as long distance nonperturbative ones, such as confinement. 

In order to relate Casimir energy calculations to correlators of the stochas- 
tic model, we will use a framework where the boundary conditions are im- 
posed by means of auxiliary fields living on the bag boundary. In this way, 
we are able to show that the lattice-adjusted SVM parameters imply that 
the simple model that computes Casimir energies treating gluon modes as 
massless photon-like gauge fields is not reliable; rather, a scenario with a 
suppressed Casimir energy due to gapped nonperturbative effective modes 
is favoured. Thus, a physical picture arises where the suppression of the 
gluon Casimir energy leaves room for the effect of center of mass quantum 
fluctuations to dominate the R~ x contribution to the bag energy. 

The structure of this paper is as follows: in section [2], as a preparation for 
the nonabelian case, we use the approach of p[5] to construct the path inte- 
gral for the Casimir energy for a confined Abelian gauge field with magnetic 
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boundary conditions, in terms of the field-strength correlation function. In 
section [HI the same approach is followed for a nonabelian gauge field within 
the context of the SVM. In section 0J we show that the same vacuum energy 
that one obtains from the SVM may be alternatively understood as due to a 
pair of fields, one of them vectorial and the other pseudovectorial, both en- 
dowed with non-standard kinetic terms. This parametrization is then applied 
to infer the main properties of the SVM Casimir energy. 
In section El we present our conclusions. 

2 Confined Abelian field 

Let us first consider an Abelian gauge field in 3 + 1 dimensions, confined 
to a spatial region V whose boundary is a regular (static) surface E. On that 
surface, we shall consider boundary conditions corresponding to a perfect 
conductor of magnetic charges. 

We shall assume that E can be described in parametric form, 

r:U — ► M 3 
(C\C 2 ) — r(C\C 2 ), (1) 

where U is a finite region of IR 2 , the domain of the parameters (}, ( 2 . In 
order the impose the boundary conditions, it is convenient to introduce first 
some notation and conventions. Since the boundary conditions are static, 
they shall hold on the world-volume r = IR x E swept by E in the course 
of its (trivial) time evolution. This world-volume can be parametrized by 
a function r M (// = 0,1,2,3), depending on C 1 , ( 2 as well as on an extra 
parameter (°: 

r^-.RxU — ► IR 4 

(C°,C\C 2 ) — r,(C), (2) 

where r (C) = Co and n(() = r f (Ci, Ca)- 

Regarding the 'target' space, we shall work in Euclidean spacetime, with 
coordinates x M , fi — 0, 1,2,3, and a metric which is the identity matrix (no 
special meaning will be given to the position -up or down- of space and 
spacetime indices). Letters from the beginning of the Greek alphabet run 
from to 2, those from the middle do so from to 3, while Roman ones 
can take values from 1 to 3. Besides, a, b, c, . . . will be reserved for surface 
coordinates, and therefore may assume the values 1 or 2. 
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At any point of T, we have the three tangent vectors ££J = which are 
not necessarily orthogonal but certainly linearly independent vectors 0. Of 
course, t\ and t 2 are tangent to E, at all times. 

Besides, we note that each component t M of t transforms as a covariant 
vector under reparametrizations: t' a ((') = ^r^(C)- 

Finally, we shall also make use of n, the (outer) normal to S. At any 
point r(C) it can be easily obtained as the vector product of the two tangent 
vectors ti and t 2 : n = t 1 x t 2 . This vector may be easily normalized by 
taking into account the fact that: 

n 2 (C) = (ti(C) x t 2 (C)) 2 = gdC) , (3) 
where g-% = det(g a b), and g a b is the induced metric on E: 

g ab (0 = (C)tf(C) ; (4) 

thus, the unit normal n is simply: n = A = -== n. As a matter of fact, 

' y J |n| ' 

since the time evolution of the boundary is trivial, gr, the determinant of the 
induced metric on T, coincides with gg. 

Regarding the transformation properties of n, one easily sees that it be- 
haves as a pseudo-scalar under surface reparametrizations. 

Now we introduce the magnetic boundary conditions: they amount to 
surrounding the region V with a perfect magnetic conductor, so that the 
tangential component of B at every point of E, and for every time vanishes: 

fi(C)xBWO) =0, (5) 

or, in terms of F^, 

^iMOK'(C) = 0, 2 = 1,2,3. (6) 

Note that they do not amount to three independent conditions, since they 
satisfy the condition: hiFijhj = 0. Thus, © may be written equivalently as: 

rj a = t\F i3 W = , a = 1,2, (7) 

where only the two independent constraints, r/i and r/ 2 appear (note that 
t' l = 0, thus there is no other constraint). 

Let us now consider the Euclidean path integral that result from imposing 
the previously introduced constraints into the vacuum transition amplitude 



x At least under the assumption that E is a regular surface. They may be normalized, 
but that step shall not be necessary. 
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for the Maxwell field. To that end, we shall define Z mag , the partition func- 
tion corresponding to a magnetic charge conductor, since it is more similar 
to the case we will need to deal with in the nonabelian case: 

z = \T)A\ P - s si A ] 

*-"mag — / i*^ ■ rl \mag ° 

r 2 

= M VA n%.]e- 5fl[A] , (8) 

a=l 

where S g is the gauge field action (including gauge fixing), and [VA] mag is the 
integration measure for the gauge field, assuming it satisfies the correspond- 
ing boundary conditions on S. Following [15], those conditions have been 
imposed, on the second line, by means of two S functionals. The integration 
measure for A^ becomes then unconstrained. 

It is convenient to exponentiate the 5 functionals, by means of two La- 
grange multipliers, denoted by 4> a ((). It has been shown that the expo- 
nentiated ^-functionals have to be invariant under reparametrizations of the 
surface [TH]. Thus, it is immediate to see that the proper expression is: 

2 r 

U*[ife] = / V<f>e- 8 **A (9) 

0=1 "* 

where 



s s [d>,A] = -i d 3 cVoA<)r(Ova(0 • (io) 



Reparametrization invariance of Ss is easily verified; indeed, d 3 ( ^/ gr(() 
is invariant by itself, while r) a is a covariant vector, whose change is com- 
pensated by the variation of the contravariant vector (fi a . Note that the 
y/gr(j£) = v^sIO factor may be cancelled with a like one in the unit nor- 
mal; thus: 



Ss[<f>,A] = -i J rf 3 C0 a (CK(C)^(r M (C))^(C) 

-i fd 3 ( (p a (C)d a rKC)Fa (r M (C))^^ fe (C) V(C) ■ (H) 



This term may be regarded as an interaction between the gauge field and a 
'current' J E due to the boundary conditions: 

S s = ~\ I d 4 xj£(x)F kl (x) (12) 
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with 



J%(x) 



d 3 ( 4> a (0 5 {4) {x„ - r M (C)) [{d a r k )m - (d a ri)n k ] ({) . (13) 



Thus, in terms of Z(J pu ), the generating functional of field strength cor- 
relators: 

p—vvyJ^v) — J_ 



Z(J f 



fVA 



>Sg(A) 



Z mag may be written as: 



i nag 



V<\> exp[-W(J^)] 



(14) 



(15) 



where J^ v agrees with (fTBl when both indices are spatial, and is assumed to 
vanish when one of them equals zero. Of course, the integral is Gaussian, 
and one is able to write the result: 



W(J ; 



\1,V J 



d A x / d A y J llv {x){F pv (x)F p(7 {y))J p(7 (y) , 



where: 



with 



and 



(F pu (x)F pa (y)) = A^ pa (x-y) 
d 4 k 



^pu,pa (•£ 



(27T) 



ik-x ^ 



pi/, pa 



(k) 



Apu,pa{k) 



1 



Ova fcufcp H~ o up hi/k a 6 no- ki/hp Ojjn knh c 



J^2 \^ uu P P PP ,v v ,v a w fia ivi/ivp ^vp ,v fi ,v a 

The coordinate-space version of (|T9l) may be expressed as follows: 
1 



*pv,pa 



where 



2 1 d 1 1 

z = x - y , D (z ) = 



DJz 2 ) 



7T 



dz 2 



7l 2 Z 4 ' 



(16) 

(17) 
(18) 

(19) 

(20) 
(21) 



The vacuum energy E in the presence of the boundary E is then given by: 

E Q = - lim ^\iaZ mag , (22) 
where T denotes the extension of the Euclidean time coordinate: \xq\ < T/2. 
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3 Confined nonabelian gauge field 

For a region surrounded by a perfect conductor of chromomagnetic charges, 
we have at the boundary the condition: 

B x n = (23) 

where B = B U T U and T u are the generators of the Lie group. The index u 
runs from 1 to N, the dimension of the adjoint representation (i.e., iV = 8, 
for 577(3)). In terms of the nonabelian field strength, 

^•MO)^-(C) = 0. (24) 

As the constraints (|24"1) are homogeneous, they can also be written in the 
equivalent form 

i^(y,x;C)%(C) = 0, VxeS, (25) 

obtained by covariantly shifting the field strength from an arbitrary 
(fixed) reference point y. F^ v (y,x;C) is the shifted field strength, defined 
by, 

F^(y,x;C) = V(y, x; C) F^(x) V~\y, x; C), (26) 
where V is the parallel transporter, from x to y, along a path C, 

V(y, x- C) = V exp f-ig J dzPA^z)) , (27) 

and V is the path ordering operator. 

The Casimir energy obtained by imposing either the constraints (l24"l) or 
the equivalent ones f l25l) should be the same, since they differ by a gauge 
transformation @. 

Thus we will impose the constraints, 

V u a = 1&F%{x,y;C) = 0. (28) 

Introducing, as in the previous section, Lagrange multipliers 0*, we obtain 
the representation, 

Z mag = J [V<f>] e~ S ^ (29) 

defined as in the Abelian case, but with 0, being a pair of fields in the adjoint 
representation, and with a generating functional that now corresponds to the 
nonabelian case: 

W{J) = -lnZ(J) , (30) 

2 This could be verified, for example, by means of an appropriate transformation of the 
corresponding Lagrange multipliers. 
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with 

Z{J) = J [Di]e^» (i4) + ^ A ^ fcx;C)jrw , (31) 

where S g (A) is the gauge-fixed Yang-Mills action. 

The vacuum energy E Q , defined as in the Abelian case, must be gauge 
invariant. This fact, naturally expected from the gauge invariance of the 
boundary conditions at S, may be verified explicitly as follows: under a 
gauge transformation associated with the group transformation U(x), the 
shifted field strength transforms in the adjoint representation, with the in- 
duced group rotation evaluated at y. This change can be compensated by a 
parallel transport (also in the adjoint representation) of the Lagrange mul- 
tipliers to the point y. Being that a unitary transformation, no non-trivial 
Jacobian is generated. 

We also note that, had an Abelian theory been considered, we would 
have F^^y, x; C) = F^ u (x) and this calculation would have reduced to the 
one performed in the previous section. 

Returning to the general nonabelian case, taking functional derivatives of 
Z( J) at J = 0, we may generate field strength correlators, 

(F£ ui (y, x x - d) . . . F% Vn (y, x n ; C n )) (32) 

which, due to colour conservation satisfy: 

(F^ vl (y, Xl ;Ci)) = 0, (33) 

and 

^{F;i{ y ,x-C)F%{y,x'-C')) = -^F^^x^y-CC) . (34) 

It can also be shown that the correlator in (1341 may only depend on C — C, 
the curve connecting x to x' . 

To proceed, we recall that the SVM is based on the strong assumption that 
this correlator does not depend at all on C and C. Moreover, the expectation 
value of the product of any number of shifted F's vanishes, while the product 
of an even number can be decomposed in terms of the sum of the different 
products that can be formed with just two field correlators. In other words, 
the assumptions of the SVM are tantamount to saying that the quantum 
fluctuations are controlled by an essentially Gaussian measure. To define the 
model, it is then sufficient to specify a quadratic W(J) functional, with the 
structure: 
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As a consequence, 

JV 



7 

■^mag 



T)(j) e f d^xd x' jJ^ v (x)F lllJ;pcr (x,x , )J^ cr (x') 



(36) 



where now refers to just a single pair of fields <fi a (without any colour index). 

The correlator for two shifted field strengths is frequently parametrized 
as follows H: 



F^.^x.x') = {5 PP 5 V(J - 5 pa 5 up )nD(-z 2 

+ - [duiza^p - z p 5 pu ) + d p {z p 5 va - z a 5 vp )\ (1 - k)£>i(-z 2 ) 

(37) 



where z = x — x' . 

Then, within the SVM, the Casimir energy is just N times the one com- 
puted in some kind of generalized Abelian theory, namely, the one which 
comes from the modified kernel F pu;pa . Note, in particular, that it contains 
a D-term, an object that is not present in the Abelian case (because of the 
absence of magnetic monopoles). This new term comes from the nonabelian 
character of the theory, and can be associated with the dual superconductor 
scenario [16] of the QCD vacuum, representing the effects of a chromomag- 
netic monopole condensate; it is crucial to derive the confining linear poten- 
tial between heavy charges. The second term (Abelian part) contains the 
modified nonlocal kernel D\, where perturbative effects dominate. We will 
discuss later specific fittings of D and D 1 which are obtained from lattice 
data. 



4 Field representation of the Casimir energy 

An equivalent functional integral representation for the Casimir energy, pre- 
viously expressed in terms of shifted field strength correlators, can be con- 
structed in terms of Gaussian fields, albeit with non standard propagators. 
We have found that the most economical parametrization is in terms of a 
vector field A p and a pseudo- vector field 

F pv = d p A u - d u A p + e pupa d p <& a , (38) 

where the new ^-dependent term represents the monopole sector of non- 
perturbative Yang-Mills. Similar field representations have been considered 

3 For a recent review see, for example, ref. [T3] . 
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to discuss Wilson loops, also in the context of the method of field strength 
correlators [17J. 

Then the (FF) correlation function is now given by: 

(F pu (x)F pa (y)) = 

= ((d^Au - d 1/ A lx )(d p A (7 - d a A p )} + e^ M v W^m'^/X^'^)^'^)) 
= O pU)P<J (G S) + Ifj,u,p<jd S , 

(39) 

where (A p (x)<S> u (y)} = 0, 

(A p (x)A u (y)) = 5,Mz 2 ) , (%(x)Mv)) = ^uS(z 2 ) . (40) 
The prime represents derivative with respect to z 2 , and we used the notation: 



1 pv , per 



o 



pv , pa 



2{d u [5 pp z a - 5 pa z p ] + d^S^Zp - S up z a ]} . 



(41) 
(42) 



Then, to reproduce the lattice parametrization of (1371) . one can make the 
identifications, 



d 2 S = kD , G' - S' = (1 - k)Z>i/4 , (43) 
which can be solved by imposing the following equations, 

S = -«^ , G=-k^ + {\-k)H , H' = D 1 /A. (44) 
In other words, the (only) nontrivial propagators are, 



(M*)Mv)) 



d A k 

(2^ 



k£>+(1- n)k 2 H 
6, 



k 2 



kD 



>VV_ ik-(x-y) 

k 2 



(45) 



We see that D and k 2 H became associated with 'structure functions' for the 
propagators of the A p and $ M fields. 

Then, the partition function (l36l) can be also represented by, Z mag = Z Aq>1 
where 

Z A ^ = J V(j)VAV^ e Se[A]-Smm + ^Jd 4 x J^(d„A v -d v A f ,+^ vptT d p ^ tT ) _ (4g) 
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Path integrating over the Lagrange multipliers 0, we get the constraints, 

Fijrf = (diAj - djAi + e ijpa d p $ a )n j = , (47) 

namely, the Yang-Mills fields Casimir energy, due to Gaussian dominance 
in the shifted field strength correlators, can be written as N times the one 
obtained for a field theory with action S e [A] + 5 m [$]. Here, both S e [A] and 
are (typically) nonlocal versions of Maxwell-like actions; for instance, 
S m [$] ~ / d^xG^D-^d^G^, with G^ = d^-d^^. We shall assume 
that a Lorentz gauge-fixing condition for A^ and $ M is implicitly included in 
the functional integral measure. 
Then, we have, 

Z A ^ = J [DAV$] mag e- s °W- s ™W ; (4 8 ) 

where the path integral is carried over fields A^, $ M satisfying the con- 
dition (j4T|) . Note that the original constraints yield boundary conditions 
corresponding to a perfect conductor of magnetic charges, while the field 
strength receives contributions from both 'electric monopole' (A) and 'mag- 
netic monopole' (<&) sectors. 



5 Analysis of the gluon contribution 
to the Casimir bag energy 

For the case of pure SU (3) @, typical ansatze for D and D\ are (see ref. [T4]): 



2\ ^2 27 _ 2 f ,4 , k 2 



24 64 

and 



D(-z 2 ) = ——or 1 \ d A k j expUkz) , (49) 

"'' 1 (P + (3vr/8a) 2 ) 4 { ' [ ' 



D^-z 2 ) = (z 2 )- 2 / dv2vD(-v) . (50) 
Jo 

We note that for large z 2 , D ~ exp ^— and, 

D.i-z 2 ) ~ (^ 2 )- 2 . (51) 



4 The effect of dynamical fermions on field strength correlators can be studied on the 
lattice, see ref. [20J. The model parameters get renormalized, however, our general analysis 
is left unchanged. 
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Although a pure SU(3) lattice calculation fixes the parameters to k — 0.74, 
a = 0.35 fm and G 2 = (496 MeV) 4 , see [TT], E] and references therein, it is 
interesting to study what happens when k goes from to 1. 
When k = 0, there is just one propagating field, 

(A,(x)A u (y)) = -H6, u , (52) 



the leading contribution to D\ is perturbative, D\ ~ l/^ 4 , and from ( l44l) 
we see that H ~ l/z 2 . Thus, the leading behaviour of the propagator 
is ~ {1/^)5^. As this corresponds to the usual behaviour of QED(4), the 
Casimir effect at k = is expected to be repulsive for a bag, and after 
multiplying by N = 8, E c « +0.7/.R. 

Then, for k close to zero, the usual result for the simplified model where 
the Casimir bag energy is computed considering the gluon modes as photon- 
like fields is recovered. Of course, the lattice fitting value n = 0.74 is far 
from that regime, and it comes as no surprise that the Casimir bag energy 
for the simplified model, after adding contributions coming from center of 
mass fluctuations, is far from the required value to adjust hadrons. 

Since k is in fact closer to one than to zero, let us analyze the opposite 
regime, when k = 1. The field propagators are now given by, 

{M*)Mv)) = (*MMv)) = -D 5j gr ■ (53) 

We can also write, 

S e [A] + S m m = \{A, (J ( 54 ) 

(where in D we have replaced k 2 by — d 2 ). 
Now, we note that (1491 implies, 



^oc k \ , A = (3vr/8a), (55) 

(Ar + A 2 ) 

so that the kinetic operator in ([54]) is proportional to, 

VW-d 2 + A 2 ) 4 , (56) 

and the Casimir energy implied by ( |48l ) is four times the one for a kinetic 
operator, 

<WW-d 2 + A 2 ). (57) 

We see that the A and $ propagators are both infrared suppressed, since 
the associated modes are in fact massive. The lattice adjusted value of a = 
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0.35 fm, corresponds to A ~ 660 Mev^ 1 . Casimir energies associated with 
massive modes at this scale, have been precisely discussed in refs. [HI fT9], 
and are strongly suppressed, for a typical hadron radius, when compared with 
contributions of order ~ 1/R. Typical values of R for the smaller hadrons 
are ~ 2/3 fm, ~ 1 fm (1 fm « ^Mev~ r ). 

It is interesting to note that the field representation we have obtained 
here is quite similar to the simple infrared modified bag model introduced in 
ref. [H]. In that reference, information associated with a mass scale A was en- 
coded in a quadratic Abelian field theory, containing an effective gluon prop- 
agator, thus discussing possible effects of nonperturbative infrared physics 
on Casimir energy calculations. 

With regard to the equivalent field theory discussed in the present work, 
three important remarks are in order: i) the field theory has been derived 
rather than assumed: its quadratic nature is a direct consequence of the 
Gaussian dominance in shifted field strength correlators, the 'hard' part of 
the calculation, namely, the generally nonlocal kernels in the kinetic terms, is 
borrowed from lattice data; ii) as it is based on field strength correlators, the 
formulation is SU(3) gauge invariant, as well as the defined Casimir energy 
and the nonperturbative information encoded in the kernels; iii) It contains a 
vector and a pseudovector field, representing electric and magnetic monopole 
sectors; the field strength receives contributions from both. 

For the lattice value of k = 0.74, the Casimir effect will be presumably 
closer to what happens at k — 1, where the effect is suppressed with respect 
to contributions of order ~ l/R, than to the value of ~ +0.7 / R, at k = 0, 
implied by "photon-like" gluons. So that the R~ l term in the bag energy is 
expected to be attractive, dominated by center of mass fluctuations —Z cm /R, 

6 Conclusions 

In order to fit properly hadron spectroscopy, the MIT bag model requires 
a term in the bag energy of the form —Z/R, with Z of order one. As it 
is well known, when Casimir energy calculations are done by means of the 
crudest simplification, considering 'photon-like' gluons, the effect turns out 
to be ~ +0.7/ R, with the wrong sign. Improvement on this disappointing 
result has been, for many years, blocked by the fact that a full nonabelian 
calculation would be extremely difficult to perform. 

In this article, we have introduced a middle way to approach the problem, 
namely, to take the most relevant nonabelian effects into account. 

Using a technique based on auxiliary fields living on the bag boundary, 
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to impose the bag constraints, we have been able to express the Casimir bag 
energy in terms of shifted field strength correlators in 577(3). Of course, 
the obtained Casimir energy representation is invariant under 577(3) gauge 
transformations. 

As lattice simulations show, the shifted field strength correlators pos- 
sess Gaussian dominance and can be parametrized in terms of two structure 
functions fitted with the lattice data. This is the so called stochastic vacuum 
model of QCD and it consistently describes short distance perturbative as- 
pects, such as asymptotic freedom, as well as large distance nonperturbative 
ones, such as confinement. Therefore, it is a natural approximation scheme 
to discuss whether the simple model that treats gluon modes as free Abelian 
massless fields can be applied, and if not, what are the consequences of taking 
nonperturbative effects into account. 

Gaussian dominance enables the introduction of an equivalent quadratic 
field representation in terms of vectorial and pseudovectorial effective fields, 
representing electric and magnetic monopole degrees of freedom, respectively. 
This representation contains two generally nonlocal kernels, whose typical 
ansatze contains a k parameter weighting the nonabelian nature of gluon 
fields. For n close to zero, the effective model is dominated by perturbative 
physics, so that in this case a "photon-like" gluon model would be obtained, 
together with the associated Casimir energy, ~ +0.7 / R. 

However, the lattice adjusted value is k = 0.74, closer to k — 1, a value 
where the equivalent Gaussian fields have a mass gap A = 660 Mev and the 
Casimir energy for typical hadron radius is strongly suppressed. Therefore, 
the Casimir energy of pure SU (3) theory is expected to be suppressed, when 
compared with attractive center of mass contributions, —Z cm /R, Z cm of order 
one. An important point is that the infrared suppression of propagators, 
in our equivalent field representation, is SU(3) gauge invariant: the fields 
represent effective modes with a nonperturbative gap, associated with SU(3) 
gauge invariant structure functions of shifted field strength correlators. 

The predicted bag energy term is, therefore, —Z/R, Z rs - / Z cm r * u X. To the 
best of our knowledge, this is the first controlled estimation of the 1/R term 
in the bag energy, and it turns out to comply with hadron phenomenology. 
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